The use of relativistic action in strong-field nonlinear photoionization 
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Nonlinear relativistic ionization phenomena induced by a strong linearly polarized laser field 
are considered. The starting point is the classical relativistic action for a free electron moving in 
the electromagnetic field created by a strong laser beam. This action has been used to calculate 
semiclassical transition rates. Simple analytical expressions for the ionization rate, the photoelectron 
emission velocity and for the drift momentum distribution of the photoelectron have been found. 
The analytical formulas apply to nonrelativistic bound systems as well as to initial states with an 

^\ • energy corresponding to the upper boundary of the lower continuum and to the tunnel as well as the 

multiphoton regime. In the case of a nonrelativistic bound system we recover the Keldysh formula 

0^ , for the ionization rate. Relativistic effects in the initial state lead to a weak enhancement of the 

rate of sub-barrier ionization and to the appearance of a nonzero photoelectron leaving velocity. 
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I. INTRODUCTION 



Relativistic ionization phenomena induced by strong laser light has become a topic of current interest |1H7| . In the 

nonrelativistic theory it is assumed that the electron velocity in the initial bound state as well as in the final state is 

•^ ■ small compared with the speed of light. However, the electrons may be accelerated up to relativistic velocities in an 

intense electromagnetic field produced by modern laser devices. If the ponderomotive energy of the electron is of the 

order of the rest energy a relativistic consideration is required. Relativistic effects in the final states become important 

f> f for an infrared laser at intensities of some 10 16 W cm~ 2 . The minimal intensity required for relativistic effects increases 

by two orders of magnitude for wavelength corresponding to visible light. The main relativistic effects in the final 

I ' state are [Q,||-0 : (i) the relativistic energy distribution and (ii) the shift of the angular distribution of the emitted 

£> , electrons towards the direction of propagation of incident laser beam. Relativistic effects have also to be taken into 

account if the binding energy Ei, in the initial state is comparable with the electron rest energy g. A relativistic 

formulation is necessary for the ionization of heavy atoms or singly or multiply charged ions from the inner shells. 

This paper is aimed to consider the relativistic effects connected with relativistic final states velocities and/or low 
lying initial states from a unique point of view. 

Let us start with the classical relativistic action for an electron of charge e moving in the field of an electromagnetic 
plane wave with the vector potential A(t — x/c). Here and below A denotes a two-dimensional vector in the y-z plane. 
The action may be found as a solution of the Hamilton- Jacobi equation and reads || 

S^ ) = m Af- r --a X -- l + a * + f \ t-^) + ^ffAdZ--^fA*dA, (1) 
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where a and / = (oi, 0,2) are constants, r = (y, z); further is £ = t — x/c, £0 is the initial value. By applying the usual 
Hamilton- Jacobi method we take the derivative of the action Sf with respect to the constants a\ , 0,2 and a and set 
the result equal to new constants /3i, P2 and 0$ in order to obtain the electron trajectory under the influence of the 
wave field. Assuming a harmonic plane wave of linear polarization with the electric field E = Fe y cosu;£ we obtain 



C^ ■ that the electron motion in the field and in the laboratory coordinate system is given by 
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where fa, fa and fa together with <X\,a% and a have to be determined from the initial conditions for position and 
velocity. Further we have introduced the notations fa 2 = 1 + a\ + a 2 , + e 2 /2, S 2 = 1 + a\ + a 2 and the parameter 
e = eF/ujmc characterizing the strength of relativistic effects. 

Consider now the process of nonlinear ionization of a strongly bound electron with a binding energy Eb comparable 
with the rest energy. Recently the ionization process in static crossed electric and magnetic fields has been considered 
[g,H . The results of this paper may be applied to the ionization in laser fields only for the case of very strong fields 
e>l, With an increasing frequency of the laser light (especially for a tentative x-ray laserWery high laser intensities 
are required to satisfy this condition. Therefore it is necessary to generalize the result of |2[g] to the case of nonzero 
frequencies. We consider the sub-barrier ionization. The condition to be satisfied is the opposite to the case of pure 
classical ionization, F <C Fb, in addition we have the quasiclassical condition Tilo <C £*. No restrictions are applied 
to the parameter e. Thus we will cover both the regime of relativistic tunnel and multiphoton ionization. 

We employ the relativistic version of the Landau-Dykhne formula MM . The ionization probability in quasiclassical 
approximation and with exponential accuracy reads 

PFcxexpj-|lm (5/(0;to) + Si(*o))} , (3) 

where Si = Eyta is the initial part of the action, Sf is given by Eq. (|l|). The complex initial time to has to be 
determined from the classical turning point in the complex half-plane 
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E f (t ) = mc 2 { i —fA(t ) + - A 2 (t ) ) = E Q = mc 2 - E b . (4) 



The minimization of the imaginary part of the action leads to the following boundary conditions |TT| 

(x,r)(t o )=0, lm(a:,r)(t = 0)=0. (5) 

In order to obtain simple analytical results we consider the case of linearly polarized laser light. Then by minimizing 
the action we obtain from Eqs. (0) and (ft) that / = 0. Further we obtain a system of nonlinear equations for the 
determination of complex initial time to and constant a, 



to = itq = arsinh ( ny 1 + a 2 — 2aeo 
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with the dimensionless initial energy e = Eq/itic 2 and the relativistic adiabatic parameter 77 = e _1 = uimcjeF. 
Substituting the values to and a into the final state action we obtain the probability of relativistic quasiclassical 
ionization in the field of linearly polarized laser light. Within exponential accuracy we get 
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where a has to be taken as the solution of Eqs. (@). Further 7 = \/2mEtuj / eF is the common adiabatic Keldysh 
parameter from nonrelativistic theory M and 7(a) = r/^/1 + a 2 — 2aeo is an a-depending adiabatic parameter. Equa- 
tion (pi) is the most general expression for the relativistic ionization rate in the quasiclassical regime and for field 
strength smaller than the above-barrier threshold. It describes both the tunnel as well as the multiphoton ionization. 
It is the relativistic generalization of the famous Keldysh result |12] . 

Consider now some limiting cases. In the limit of tunnel ionization 77 <C 1 we reproduce the static result of Refs. 
P,p| and obtain the first frequency correction 

W oc exp I -$ 

#= 2V3(l- a yA_3V3(l- a ^ 
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where F s = m 2 c 3 /eh = 1.32 • 10 16 V/cm is the Schwinger field of quantum electrodynamics E3| and ao = (eo + 
y/eQ + 8)/4. In the nonrelativistic regime, £& = Eb/mc 2 <§; 1, the parameter ao = 1 — £f,/3 + e 2 /27 and the probability 
of nonrelativistic tunnel ionization including the first relativistic and frequency corrections reads 



W ex exp < — - 
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Here the first two terms in the brackets describe the familiar nonrelativistic ionization rate including the first frequency 
correction B, the next two terms are the first relativistic corrections. It follows from Eq.(|8j) that the account of 
relativistic effects increases the ionization rate in comparison with the nonrelativistic rate. However, even for binding 
energies of the order of the electron rest energy the relativistic correction in the exponent is quite small. In the 
"vacuum" limit Eq. (0) results into W ex exp{— 9Fs/2F(l — 9/40?7 2 )}. We find a maximal deviation of about 18% 
in the argument of the exponential from the Keldysh formula. Here the "vacuum" limit shall not be confused with 
the pair creation from the vacuum. It is known that there a no nonlinear vacuum phenomena for a plane wave [pL3|| . 
In contrast to that we deal here with the ionization of an atom being in rest in the laboratory system of coordinates. 
We also mention that we employ the single particle picture. Therefore the pair production processes are beyond the 
scope of the present paper. 

Consider now the multiphoton limit 77 ^> 1 . In this case the parameter a = 1 — £b/21n27 and the ionization 
probability in the relativistic multiphoton limit reads 
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Again the first two terms in the brackets reflect the nonrelativistic result |0| , the relativistic effects which lead to an 
enhancement of the ionization probability are condensed in the third term. 

It has been shown that there is an enhancement of ionization rate in the relativistic theory for both large and small 
77. This should be compared with the results found by Crawford and Reiss. In their numerical calculations they also 
found an enhancement of relativistic ionization rate for a circularly polarized field and for 77 3> 1, but for 77 <C 1 
their results suggest a strong reduction of the ionization probability H . For the case of linearly polarized light the 
ionization rate is found to be reduced by relativistic effects ||. However, Crawford and Reiss studied the above-barrier 
ionization of hydrogen atom within the strong-field approximation. In contrast to that we have investigated the sub- 
barrier ionization from a strongly bound electron level, which yields an enhancement of the ionization rate. This 
enhancement is connected with a smaller initial time to- As a result the under barrier complex trajectory becomes 
shorter and the ionization rate increases in comparison with the nonrelativistic theory. Figure 1 shows the relativistic 
ionization rate Eq. (R) and the nonrelativistic Keldysh formula as a function of the binding energy e&. The figure 
should be considered only as an illustration of the enhancement effect. The frequency and intensity parameters used 
for the calculations are still not available for the experimentalists. 

The switch from the multiphoton to the tunnel regime with increasing field strength may be studied in the nonrel- 
ativistic limit Eb <C 1. Here within first order of £&, with a = 1 — (£{,/27 2 )[(7/arsinIi7)-y/l + 7 2 — 1], the ionization 
probability is found to be 
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The terms in /(y) which do not vanish as Eb — > represent the nonrelativistic quasiclassical ionization rate foun d by 
Keldysh |1J] ; the terms proportional to £& are the first relativistic correction to the Keldysh formula. Equation (|ll| ) 
is valid in the whole 7-domain, i.e., in the multiphoton regime 7 < 1 as well as in the tunnel limit 7 > 1. For small 
adiabatic parameters, i.e., 7 — ► 0, it coincides with Eq. (g); in the case of large 7 — ► 00 it transforms to Eq.(|Hj). We 
mention that Eq. (TO]) reproduces the full relativistic formula Eq. (R) with very high accuracy for Eb < mc 2 . 

Consider now the modifications of the energy spectrum induced by relativistic effects. In the nonrelativistic theory 
and in the case of linear polarization the most probable value for the electron momentum at the time of emission, 
t = 0, is zero. The electrons are preferably emitted in the direction of the polarization of the laser beam. In the 
relativistic theory employed in this paper we may set the constants a\ = 02 = in Eqs. (0). Then we obtain for the 
most probable emission velocity in the laboratory system of coordinates 

1 - a 2 
v x = c— — - , Vy = v z = , (12) 

1 + or 

where a has to be taken as the solution of the second equation of Eqs. (pi). In the static limit w->0we reproduce the 
results of Mur et al. H . It follows from these equations that a strongly bounded electron is emitted in the direction of 



the laser beam propagation, i.e., perpendicular to the direction of the laser beam polarization. For a nonrelativistic 
initial state, ej < 1, the mean emission velocity along the beam propagation v x = ce(,/3 is small. Nevertheless, 
the mean emission velocity seems to be the most sensitive measure of the appearance of relativistic effects in the 
initial states. In Fig. 2 the ^-component of the leaving velocity is plotted versus the binding energy of the initial 
state. Though we have choosen the same parameters of the laser beam as in Fig. 1 it should be mentioned that 
the dependence of the emission velocity x-componcnt on the laser parameters is rather weak. The main parameter 
determining the leaving velocity along the propagation of the laser beam is the binding energy of the atom. 

The electron energy spectrum is also influenced by relativistic final states effects. We put a\ — py^/mc, <X2 = p z .o/ m c 

and a = (—p x ,o + \ 1 + Px o +-P« o +-P2 o)/ mc - The calculations will be restricted to the tunnel regime 7 <C 1. 
Assuming weak relativistic effects in the initial and final states, £5 <gc 1 and p y> a, Pz.o "C mc one obtains 
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where W is the total ionization rate Eq. (|9|) in the weak relativistic tunnel regime. The first exponent in Eq.(|13|) 
describes the momentum distribution in the plane perpendicular to the polarization axis. There is only one relativistic 
effect in the weak relativistic regime considered here - the appearance of the mean momentum at the emission time 
< Px.o >= Eb/3c. The nonzero mean emission velocity along the propagation vector destroys the symmetry in the 
(x,z)-plane that exists in non-relativistic theory. The first term in the second exponent of Eq. ( |l3l) determines the 
nonrelativistic energy spectrum for the low energetic electrons moving along the polarization axis p? < 47 2 m 2 c 2 /3, 
whereas the second, relativistic term becomes important for the high energy tail p 2 > 47 2 m 2 c 2 /3. It is only in the 
case of small adiabatic parameter 7 < 0.1, that the high energy condition does not contradict the condition jvo < vac. 
We mention that the second term in the second exponent agrees with a corresponding term of Krainov (qM . 

In conclusion, the expressions obtained in this paper within exponential accuracy may be improved by taking into 
account the Coulomb interaction through the perturbation theory. The results of this paper may be also used in 
nuclear physics and quantum chromodynamics. 
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FIGURE CAPTIONS 

(Figure 1) Absolute value of the logarithm of the ionization rate — In W versus the binding energy of initial level 
e& = E^/mc 2 . The solid line shows the relativistic rate Eq.(Q), the dashed line is the nonrelativistic Keldysh 
formula (Eq. (11) without the relativistic correction term). The curves are shown for a frequency to — 100 and 
an intensity / = 8.5 • 10 7 (in a.u.). 

(Figure 2) The ^-component of the emission velocity v x /c versus the binding energy of initial level e?, = E^/mc 2 . 
The emission velocity in the nonrelativistic theory is zero. The curve is shown for a frequency w = 100 and an 
intensity / = 8.5 • 10 7 (in a.u.). 
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